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Abstract 

The dielectric susceptibility of a wide class of dielectric materials follows, over ex- 
tended frequency ranges, a fractional power-law frequency dependence that is called the 
"universal" response. The electromagnetic fields in such dielectric media are described 
by fractional differential equations with time derivatives of non-integer order. An exact 
solution of the fractional equations for a magnetic field is derived. The electromagnetic 
fields in the dielectric materials demonstrate fractional damping. The typical features 
of "universal" electromagnetic waves in dielectric are common to a wide class of ma- 
terials, regardless of the type of physical structure, chemical composition, or of the 
nature of the polarizing species, whether dipoles, electrons or ions. 
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1 Introduction 



A growing number of dielectric relaxation data show that the classical Debye behavior pQ is 
hardly ever observed experimentally [21 El H]. The fact that different dielectric spectra are 
described by the power laws is confirmed in dielectric measurements realized by Jonscher 
[21 E] for a wide class of various substances. The dielectric susceptibility of most materials 
follows, over extended frequency ranges, a fractional power-law frequency dependence, which 
is called the law of "universal" response [21 [3]. This law is found both in dipolar materials 
beyond their loss-peak frequency, and in materials where the polarization arises from move- 
ments of either ionic or electronic hopping charge carriers. It has been found [51 [6] that 
the frequency dependence of the dielectric susceptibility x(u) = x'( w ) ~~ ^x"^) follows a 
common universal pattern for virtually all kinds of materials. Namely, the behavior 

X '{u)~u n -\ (0<n<l, u;»u; p ), (1) 

and 

X '(0) - x'M ~ ""N X"(u)~u m , (0<m<l, u^u p ), (2) 

where * s the static polarization and uo p the loss-peak frequency, is observed over many 
decades of frequency. Note that the ratio of the imaginary to the real component of the 
susceptibility is independent of frequency. The frequency dependence given by equation ([I]) 
implies that the real and imaginary components of the complex susceptibility obey at high 
frequencies the relation 

- cot I — J , (lu » Up). (3) 



X '(u) V 2 J 

The experimental behavior (T5]) leads to a similar frequency-independent rule for the low- 
frequency polarization decrement: 

_|lJ_ =tM1 (_), (,«. p) . (4) 

There are many models that describe and explain the universal response laws. Let us note 
some of them. A general equation for the susceptibility of disordered systems is proposed in 
[7]. It is shown [7] that this equation is a good frequency domain representation of the time 
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domain Kohlrausch- Williams- Watts stretched exponential. In [8], the authors have derived 
a general two-power-law relaxation function for heterogeneous materials using the maximum 
entropy principle for nonextensive systems. The asymptotic power-law behaviors coincide 
with those of the Weron generalized dielectric function derived in the stochastic theory from 
an extension of the Le'vy central limit theorem. These results [B] are in agreement with the 
Jonscher universality principle. The memory function approach and scaling relationships can 
be used [9] as a basis for the dynamic model of symmetric dielectric spectrum broadening. 
For the correspondence between the relaxation time, the geometrical properties, the self- 
diffusion coefficient and the exponent of power-law wings was established in [9]. In [10], a 
memory function equation and scaling relationships were used for the physical interpretation 
of the Cole-Cole exponent. The correspondence between the relaxation time, the geometrical 
properties, the self-diffusion coefficient and the Cole-Cole exponent was considered in [10J. 
The fractional power laws frequently observed in dielectric measurements can be interpreted 
[TT] in terms of regular singular points of the underlying rate equation for the process giving 
rise to the dielectric response. The review [3] presents a wide-ranging broad-brush picture 
of dielectric relaxation in solids, making use of the existence of a "universality" of dielectric 
response regardless of a wide diversity of materials and structures, with dipolar as well as 
charge- carrier polarization. 

The theory of integrals and derivatives of non-integer order goes back to Leibniz, Liouville, 
Riemann, Grunwald, and Letnikov [T2J [13] . Fractional analysis has found many applications 
in recent studies in mechanics and physics. The interest in fractional equations has been 
growing continually during the last few years because of numerous applications. In a short 
period of time the list of applications has become long (see for example p3 1fT5[fT6llT71[T8lfT9] ). 
Note also the application of fractional calculus to the problems of classical electrodynamics 
[201 EQ. 

In [22] fractional calculus has been undertaken in order to understand the nature of one 
type of nonexponential relaxation. An equation containing operators of fractional integration 
and differentiation is obtained and solved [22], which the relaxation function obeys in this 
case. In [23], the potential of the fractional derivative technique is demonstrated using the 
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example of the derivation of all three known patterns of anomalous, nonexponential dielectric 
relaxation of an inhomogeneous medium in the time domain. It is explicitly assumed [23] 
that the fractional derivative is related to the dimensionality of a temporal fractal ensemble 
in a sense that the relaxation times are distributed over a self-similar fractal system. In [24J , 
the model of dielectric relaxation based on the fractional kinetics containing the complex 
power-law exponents was used to received a confirmation in description of the real part of the 
complex conductivity by the fitting function. The coincidence between the suggested model 
and measured data gives [24J the possibility to suggest a more reliable physical picture 
of the swelling process that takes place in neutral and charged gels. The real process of 
dielectric relaxation during the polymerization reaction can be described [25J in terms of 
the fractional kinetic equations containing complex-power-law exponents. Based on the 
physical and geometrical meanings of the fractional integral with complex exponents there 
is a possibility [25J of developing a model of dielectric relaxation based on the self-similar 
fractal character of the averaged microprocesses that take place in the mesoscale region. 

In this paper, we prove that in a time domain the fractional power-law frequency depen- 
dence gives differential equations with derivatives and integrals of non-integer order. We ob- 
tain equations that describe "universal" electromagnetic waves for such dielectric materials. 
The power laws are presented by fractional differential equations such that the electromag- 
netic fields in the materials demonstrate "universal" fractional damping. An exact solution 
of the fractional equations for magnetic field is derived. The electromagnetic fields in the 
dielectric materials demonstrate fractional damping. The suggested fractional equations are 
common (universal) to a wide class of materials, regardless of the type of physical structure, 
chemical composition or of the nature of the polarizing species, whether dipoles, electrons 
or ions. A possible link between dielectric response and structure of low-loss dielectrics is 
discussed. 
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2 Fractional equations for universal laws 

For the region uj ^> u p , the universal fractional power law ([1]) can be presented in the form 

X(u) = X a ^)- a , (0 < a < 1) (5) 
with some positive constant Xa and a = 1 — n. Here 

(ioo) a = |w| a e iair ' fln(w)/2 . 

It is easy to see that relation ([3$ is satisfied. 
The polarization density can be written as 

P(t,r)=F- 1 (P(w,r)) =e F- 1 (x(w)E(w, r)) = eoXa-F" 1 (M^E^r)) , 

where P(o>, r) is a Fourier transform JF of P(t, r). 

Note that the Fourier transform JF of the fractional Liouville integral [121 IS] 

(I°f)(t) = —f t mdt ' 
is given by the following result (see Theorem 7.1 in [12] and Theorem 2.15 in [13]): 

Using the fractional Liouville integral, the fractional power-law ((Sj) for x(u;) m the fre- 
quency domain gives 

P(t, r) = eoXa (/?E)(t, r), (0 < a < 1). (6) 

This equation means that polarization density P(t,r) for the high-frequency region is pro- 
portional to the fractional Liouville integral of the electric field. 

For the region to -C uj p , the universal fractional power law (jlj) can be presented as 

X(cu) = x(0)-X^f, (0 < /? < 1) (7) 

with some positive constants xp-, x(0)> an d (3 = m. It is not hard to prove that equation (j3J) 
is satisfied. 
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Note that the Fourier transforms JF of the fractional Liouville derivative [T2l [T3] 



are given by the following result (see Theorem 7.1 in [12] and Theorem 2.15 in [13]): 

{FDif){u,) = {iu>)' , {rf){u>). 

Using the definition of the fractional Liouville derivative, the fractional power law (J7J), 
gives the polarization density 

in the form 

P(t, r) = eox(O) E(t, r) - e X/3 P?E)(t, r), (0 < /3 < 1). (8) 

This equation means that polarization density P(t, r) for the low- frequency region is defined 
by the fractional Liouville derivative of the electric field. 

Equations (jSJ) and (jBJ) can be considered as the universal laws in the time domain. These 
equations allow us to derive fractional wave equations for electric and magnetic fields. 



3 Universal electromagnetic wave equation 

Using the Maxwell equations, it is easy to obtain the equation 

9 2 E(t,r) <9 2 P(t,r) 1 , , , „ dUt,r) 

For the region u ^> u p , the polarization density P(t, r) is related with E(£, r) by equation 
(EJ). Substituting (jHJ) into (JSJ), we obtain the fractional equation for electric field 

^^^ + ^(Dl--E)(t,r) + (graddzvE-V 2 E) = -^^^, (0 < a < 1), (10) 

where t> 2 = l/(e ^)- Note that divE 7^ for p{t,r) = 0. 
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For the region u <C u p , the fields P(t, r) and E(i, r) are connected by equation (JHJ). Then 
equation ([9]) gives 

^^-^P+ +/3 E) + (^arf^E-V 2 E) = - / ,|, (0</3<l), (11) 

where 

.2 1 „ _ Xfi 



e o /i[l + x(0)]' P 1 + X(0) 
Equations (TXOT) and f TTTT) describe the time-evolution of the electric field in the dielectric 

materials. These equations are fractional differential equations. 

Using the Maxwell equations, it is not hard to obtain the equation for a magnetic field 

d 2 B(t,r) 1 , , 1 d m . . 1 ,.. , 

= V 2 B (t,r ) + -—curlP {t,r ) + -curl} (t,r . 12 

or^ e fi e at e 

The experimental applied field B(i, r) can be presented as 



B(t,r) 



0, t<0, 
B(t,r), t > 0. 



For the region u ^> the polarization density P(t,r) is related to E(t, r) by equation (jSJ), 
and we obtain the fractional equation for magnetic field 

- f B ^ r) + ^ (oA 2 ~ a B) (t,r) - V 2 B(t,r) = p curl j(f, r), (0 < a < 1), (13) 

where v 2 = l/(e n), and o-D+~ a is the Riemann-Liouville derivative [H] on [0, oo) such that 

(^/IW^Sf^ (0<»<1). 



T(a)dt 2 J (t-t'y~ a 
For the region u <C we obtain 



\ 92B ^ r) ~ % (oA 2+/3 B) (t,r) - V 2 B(t,r) = ficur/jftr), (0 < (3 < 1), (14) 
Vp at Vp \ / 

where 

2 = 1 Xf3 

Vfi e ofi [1 + x(0)] ' 0/3 1 + X(0)' 
Equation (IT41 is a fractional differential equation that describes magnetic field in dielectric 

media. 
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4 Fractional damping of magnetic field 

Let us obtain a solution of equations (fT3l) and (fT4|) . These equations can be represented in 
a general form. This general form of the fractional equations for magnetic field is 



( a D?B)(t,r) - \ x ( a DfB) (t,r) - A 2 V 2 B(t,r) = f(t,r), (1 < (3 < a < 3), 



(15) 



where B(a, r) = 0, and the curl of current density of free charges is considered as an external 
source term: 

f(t,r) = /iA 2 cur/ j(t, r). 
Equation (fl5|) gives equation (fl3j) for a = 2, 1 < /3 < 2, and 

-Xa, A 2 = v 2 = l/(e n). 



Ai 



Equation (fH|) can be derived from (TT5T) with 2<a<3,/? = 2, and 



Ai 



1 



i + x(o) 



A 2 = - 



ap X/3 ap SqUXp 

An exact solution of ([175]) can be presented in terms of the Wright functions 
Using the three-dimensional Fourier transform of equation (TT5]) with respect to coordinates 
and Theorem 5.5 in [TBI, we obtain the solution 



B(t,r 

where 

OO / 

G(t-t',k) = J2 — 



1 



(2vr) ; 



cffc / dt' e* kr f(t',k)G(t-t',k), 
Jo 



k 2 x 2 y 



si 



(t - t') as+a ~ l 1^1 



[s + lX 



Here, i^/i is the Wright functions 



(as + /3, a) 



(as + a, a — (3) 

T(s+j + l) z* 
T(as + (3 + aj) j! 



Xi(t-t 



/\a-/3 



(16) 



(17) 



that can be derived as d s E a ^[z]/dz s of the Mittag-Leffler function E a> p[z] (see [T3l 26J). 
Equation (1161) describes the fractional field damping of magnetic field in the dielectric media. 
An important property of the evolution described by the fractional equations for electromag- 
netic waves is that the solutions have power-like tails. 
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5 Possible link between dielectric response and molec- 
ular structure of dielectric 



An interesting problem of considerable practical importance is the link between the dielec- 
tric loss of low-loss materials and the identifiable physical and structural features in the 
material. The following questions relating to this subject can be formulated. What material 
conditions must be satisfied for the low-loss property to be observed? What are the physical 
interpretations of this type of behavior? 

Jonscher has suggested in j2j [23 EQ] that the key to the understanding of low-loss be- 
havior lies in dipolar screening which leads to interactions between dipoles in materials with 
sufficiently low dipole densities. In other words, the Jonscher 's proposition is that elec- 
trostatic interactions between dipoles even at large distances from one another lock these 
dipoles so that their response to external alternating fields is weakened. The physical basis 
of low-loss dielectrics can be connected j2l ESS ED] with the effect of dipolar screening which 
produces interaction between neighbouring dipoles, thereby preventing a significant part of 
them from following external alternating fields. It is evident that the case of complete screen- 
ing must correspond to Debye relaxation, because there is no interaction between any of the 
dipoles. The result of dipole-dipole interactions with incomplete screening extending over 
many dipoles is that the system becomes insensitive in the sense that the fields generated 
by the interacting dipoles are stronger than externally applied fields [30J. 

We can also assume that the effect of screening is very important for dipolar relaxation. 
In [30], the Debye screening of the static fields is considered, but the form of dipole-dipole 
interaction is not discussed. In this section, we consider a possible form of screening dipole- 
dipole interaction for low-loss dielectric materials. In general, the dipole feels an electric 
field both from the permanent charges of the system and from the other induced dipoles. To 
describe low-loss dielectrics, we suggest using a modification of the molecular polarization 
method. This polarizable point dipoles method has been applied to a wide variety of atomic 
and molecular systems, ranging from noble gases and water to amorphous materials and 
to proteins. Let us note the main features of the polarization methods and a possible 
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generalization to describe low-loss dielectrics. In the point dipoles method, a polarizability 
a is associated to one or more sites [33]. This method for treating polarizability is to add 
point inducible dipoles on some or all sites. 

The total electric field acting on each site is produced by the external charges and by 
induced dipole moments [2B E21 E31 EH ESI HD] , such that 



E„ — E° — TnmPm, (18) 

m^n 

where p m is a dipole moment, T nm denotes the dipole field tensor, and E° is the static field 
at site n due to the permanent charges 



K = Y. ^> ( 19 ) 

ra+n ' " nm 



where q m is the charge at site m and r nm = \r nm \ is the distance between n and m. The 
induced dipoles interact through the dipole field tensor 

T nm = -^-1 - -|-r ® r, (20) 

' nm nm 

where I is the unit tensor (identity matrix). Here r <8> r is the tensor with elements XkXi, 
where Xk (k,l = 1,2,3) are the Cartesian components of the vector r nm between n and m. 
The dipole moment, p n , induced on a site n is proportional to the electric field E n at that 
site, 

Pn ^nE n . 

The proportionality constant is the polarizability tensor a n . The induced point dipole on 
site n can be computed iteratively until a given threshold of convergence for the induced 
dipole is reached (for example, see [31] for issues concerning the efficiency of the different 
methods for liquid state simulations). The energy of dipole-induced dipole interaction has 
the form 

Upp ~ ^ ^ ^ ^ Pn^-nmPm' 

m mj^n 

Here we consider that all the dipoles will interact through the dipole field tensor. The 
method of Applequist, Carl, and Fung, [33, 35J for calculating molecular polarizabilities uses 
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this approach. One problem with coupling all the dipoles with the interaction given by 
equation fl20l) is the "polarization catastrophe". As pointed out by Applequist et al [33] 
and Thole [36], the polarization p ra , and therefore the induced dipole moment, may become 
infinite at small distances. The polarization catastrophe is avoided by screening (attenuating) 
the dipole-dipole interaction [SH] . As with the screening of the static field, screening of the 
dipole-dipole interaction can be physically interpreted as correcting for the fact that the 
electronic distribution is not well represented by point charges and dipoles [33j [361 [37J . The 
Thole approach for screening is to introduce the screening (damping) functions. 

Equations ( fl8l) - (1201) retain their validity [40J with the only change being that both the 
electric field created by a fixed charge and that created by a point dipole (depending on the 
molecular model) are screening by functions fi(r) and f2(r), 



E: 



and 



1 



/l( r nm)"5 — I — /2( r nm) 



(21) 



(22) 



In the limit of point charges and point dipoles, we have /i(r) = /2O") = 1, and the usual 
expressions (fT9|) and ([20]) are obtained. If, on the contrary, they are thought to be spatially 
extended, the form of the screening depends on the charge distribution assumed. Thole 
concluded that a linear decrease of charge density (up to a cutoff a) was ideal for the 
purpose of fitting of the molecular polarizability for which it had been designed. In this 
approximation we have for the screening functions 



/i(r) = 4 



r 


3 


r 




-3 




a 




a 



/ 2 (r) 



(23) 



The possible alternative is an exponential distribution [38], [39], HO] . Note that various com- 
puter simulations with different screening functions have been used. 

We assume that Thole linear decreasing should be replaced by fractional power-law de- 
pendence for low-loss dielectrics. Therefore we suggest that dipolar screening in low-loss 
dielectric materials can be described by the screening functions 





r 


ai 


r 


«2 




— «1 






a 




a 



a 2 



Mr) 



a 2 



(24) 
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where a,\ and a.2 are non-integer positive numbers that define the fractional parameters 
a = 1 — n and f3 = m of laws ([!]) and ([S]). The simplest possible dependences are 



a\ = 3 — a, a.2 = 4 + (3. 



(25) 



For ai = 3 and a 2 = 4, equations (124"|) give the Thole screening functions fl23l) . Using 
functions (Ell) , we obtain the field from the permanent charges in the form 



E 



E 

m^n 



Q> Trim 



«i-2 



a 2 -2 



(26) 



and the induced dipoles interact through the dipole field tensor 



m^n 



r«2 


Trim 


01 


la 3 


a 





I 



air"l m l + 3r <g> r 



s~i 3 iyi2 

nm 



Trim 


0-2— 3 ^ 


a 





(27) 



For «! = 3 and «2 = 4, equations ( l26il and ( 1271) give the usual Thole screening. For treating 
polarizability of low-loss dielectric, we suggest use the fractional values of the parameters ax 
and «2- A possible generalization of exponential screening is discussed in the Appendix. 

Our suggestion to use fractional power laws in (12411 is based on the mathematical results 
that have been obtained in [41] . In this paper, we consider the vectors r n that define sites 
as in solid state physics [42], i.e., 

3 

i=i 

where aj are the translational vectors of the lattice, then 



Note that this representation can be used as an approximation, and, in general, should be 
generalized for amorphous materials. As a result, equations (1261) and (1271) describe power-law 
long-range interactions of the type \n — m| _s_1 with fractional values of s. In the continuous 
limit, these interactions are described by fractional integro-differentiation [41J. Using the 
results of [41], we can state that dipole-dipole interactions with fractional screening are 
connected with integro-differentiation of non-integer order. 
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As a result, we assume that the link between the dielectric loss of low-loss materials 
and the physical and structural features in the material can be realized by the generalized 
polarizable point dipoles method with fractional screening in molecular dynamics and Monte 
Carlo simulations. Using these simulations the presence of identified defects in the structure 
of low- loss materials can be also taken into account. To realize this suggested approach 
additional investigations are required. 

6 Conclusion 

In conclusion, the electromagnetic fields and waves in a wide class of dielectric materials 
must be described by fractional differential equations with derivatives of order 2 — a and 
2 + j3, where < a < 1 and < /3 < 1. The parameters a = 1 — n and (3 = m are 
defined by exponents n and m of the "universal" response laws for frequency dependence 
of the dielectric susceptibility. A remarkable property of the dynamics described by the 
fractional equations for electromagnetic fields is that the solutions have power-like tails. 
The typical features of "universal" electromagnetic waves in dielectric are common to a wide 
class of materials, regardless of the type of physical structure, chemical composition or of the 
nature of the polarizing species, whether dipoles, electrons or ions. We assume that the link 
between the dielectric loss of low- loss materials and the physical/structural features in the 
material can be obtained by the polarizable point dipoles method with fractional screening 
of dipole-dipole interactions in the molecular dynamics and Monte Carlo simulations. 

For small fractionality of a (or /3), it is possible to use e-expansion [27] over the small 
parameter e = a (or e = 1 — /3). There are several numerical methods to solve fractional 
equations (see for example |28j). Note that the suggested fractional differential equations, 
which describe the electromagnetic field in dielectric media with power-law response, can be 
solved numerically. For example, the Grunwald-Letnikov discretization scheme [12] can be 
used to compute fractional equations for electromagnetic field in dielectric. 
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Appendix 



In section 5, we consider the Thole screening functions and their fractional generalizations. 
The alternative distributions that are used in the polarizable point dipole method are the 
exponential screening functions 



A(r) = l-exp{-(£) 3 } (2f 



If we consider the modification of fi(r) of the form 



j 2 (r). I I exp{-(^) 3 j. (29) 



/i(r)= (l-exp<j-6('n 3 )l expj-c-k (30) 



then for r-Ca and b = 4, c = 3/4, we obtain the Thole function for the electric field created 
by fixed charge. For b = 1 and c = 0, equation (1301) gives (128]) . For r ^> a, we can use the 
approximation 

A(r) = exp{-~~} (31) 

that is the Debye screening function with = 4a/3 used by Jonscher in [21 [221 130] . 

Note that the suggested fractional screening function f\(r) of the form (124]) can be an 
approximation of the screening function 

/A)=fi-e*p{-&mw-^-r~"), (32) 



«2 v <2> 

If a\ = 3 and «2 = 4, then equation ( 1321) gives (j3~0l . For 6 = «2 and c = a\jai with 
non- integer values of «i and «2, we have the fractional screening function (124"I) . 

As a result, the exponential functions can be used in the generalized polarizable point 
dipoles method for molecular dynamics and Monte Carlo simulations of low-loss dielectric 
materials. 
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